Compaction in Granular Solid Hydrodynamics 



O 

(N 

> 
O 

On 



O 



C 

O 

i o i 

(N 
> 

(N 



On 
O 



X 



Yimin Jiang 1 and Mario Liu 2 

'Central South University, Changsha 410083, China 
2 Theoretische Physik, Universitat Tubingen, 72076 Tubingen, Germany 
(ftDated: November 22, 2010) 

Compaction is considered and embedded into broader granular behavior. Reversible compaction 
is related to the pressure exerted by agitated grains, a quantity relevant to dense flow. Irreversible 
compaction is derived from the loss of elastic deformation, the physics behind elasto-plastic flows. 
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Granular compaction - density increase by gentle per- 
turbation of the grains - is a ubiquitous phenomenon, 
relevant in many contexts such as efficient packaging, 
or when establishing the structural stability of off-shore 
windmills that are subject to continual surfs. Com- 
paction has been studied widely but is not yet fully un- 
derstood, see the review [l| , which carefully considers the 
phenomenon and cites many useful references. There are 
two compaction branches: In the irreversible one, the 
density increases monotonically; in the reversible one, 
one observes compaction and loosening, depending on 
whether the perturbation is weakened or strengthened. 

Granular solid hydrodynamics (gsh) is a set of partial 
differential equations Q constructed to model granular 
medium in all its facets, though qualitatively at present. 
It was first employed to calculate static stress distribution 
for various geometries, including sand piles, silos, and 
point load, achieving results in agreement with observa- 
tion [3] . It was then employed to consider slowly strained 
granular solid 0|, and found to yield response envelopes 
similar to those from the hypoplastic theory Re- 
cently, the critical state - generally considered a hallmark 
of granular behavior @ - was identified as a steady-state, 
elastic solution of GSH Q. Uniform dense flow, fluidiza- 
tion and jamming were also considered Q, finding broad 
agreement to experiments, and more narrowly focused 
existing theories 0, Hoj ]. Finally, the velocity of elastic 
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waves were calculated as a function of the stress 

Here, we use GSH to consider compaction. At present, 
"geometric frustration models" [l| are employed to ac- 
count for the irreversible branch, and the Edward en- 
tropy SEd 12| for the reversible one. Both are again 
narrowly focused. Moreover, the notion of SecI entails 
hidden assumptions that merit closer scrutiny. 

Derived in compliance with general principles, GSH 
consists of • five conservation laws for the density re- 
spectively of the energy w, mass p, and momentum pvt, 

• an evolution equation for the elastic strain uy, and 

• balance equations for two entropy densities, the true 
and the granular one, s and s g . Two entropies are neces- 
sary because of granular media's two-stage irreversibility: 
Macroscopic energy, kinetic and elastic, dissipates into 
mesoscopic, inter-granular degrees of freedom, mainly 



granular jiggling and the collision-induced, fluctuating 
elastic deformation. After a characteristic time, the en- 
ergy degrades further into microscopic, inner-granular 
degrees of freedom, especially phonons. The two en- 
tropies account respectively for the energy of the meso- 
and microscopic degrees of freedom. 

Now a brief presentation of the GSH-equations, tai- 
lored for the purpose of granular compaction. (See 
for the complete theory.) An important quantity of 
GSH is the energy density w, that in its local rest 
frame is a function of the above listed four variables, 
Wo{p,Uij , s gi s). The conjugate variables: true tempera- 
ture T = dwo/ds, elastic stress 7r.y = —dw /duij, chem- 
ical potential p = dwo/dp, granular temperature T g = 
dwo/dsg, and gaseous pressure Pr = pp+s g T g + sT—wo, 
are given once wq is. The Cauchy stress is given as 
(1 — a)irij +PrSij, [a is explained below Eq ((2J,] of which 
only the trace is needed here, 



P = (1 - a)7t«/3 + P T . 



(1) 



In GSH, the elastic strain u%j is taken as the portion of 
the strain that deforms the grains and leads to reversible 
storage of elastic energy. The energetically indifferent 
rest, such as rolling and sliding, is the plastic portion. 
Again, only A = — uu is relevant. Its evolution, taken as 



d t A + (1 - a)v u = -A/r, 1/r = XT g 



(2) 



accounts for two facts: First, only the portion 1 — a of 
the strain rate Vij = \i^7iVj + VjUj) deforms the grains, 
changing A and the energy u>o(A) - the rest is plastic. 
That the same a appears in Eq is a result of the 
Onsager relation. But it may also be understood via an 
analogy to a bicycle gear shift: A bigger a implies more 
pedaling at reduced torque, both with the same factor. 
Second, when grains jiggle and T g is finite, the grains 
briefly lose contact with one another. The result is a 
slow relaxation of granular deformation - as expressed 
by A/t. The full tensorial version of Eq ([2|) is closely 
related to hypoplasticity, and yields a realistic account 
of the elasto-plastic dynamics of granular media, see 

Because Eqs (|1|2|) are fairly general, it is the energy 
density w that encodes the material-specific information. 
The expression we consistently employed to study static 
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stress distribution and hypoplastic deformation is: wq 
— w%(Uij,p) +W2(s g , p)+ws(s, p), where Wi is the macro- 
scopic, u>2 the mesoscopic, and W3 the microscopic con- 
tributions. We take them (if without shear strain and 
away from the virgin consolidation line) as 

Wl (^,p) = l?A 2 ' 5 /2.5, w 2 (s g ,p) =s 2 g /(2pb), (3) 

and W3 = ep/m. Explanations: • The elastic energy w\ 
contains the lowest order term in A. Being 2.5 rather 
than two, the exponent accounts for the Hertz-like con- 
tacts between the grains and the Coulomb yield [l3| . As 
granular systems are slightly stiller at higher densities, 
the elastic coefficient B grows slowly and monotonically 
with p. In the present context, however, we may approx- 
imate B as a constant — of around 8 GPa for river sand 
and 7 GPa for glass beads Q - because the density de- 
pendence of W2 is the dominant one. • W2 contains the 
lowest order term in an expansion in s g . The exponent is 
two because W2 has a minimum at s g — 0. The density 
dependence of b(p) = b (l — p/p cp ) a , with a = 0.1, is cho- 
sen to reproduce the observed, volume-dilating pressure 
contribution from agitated grains, Pt ~ w 2 /(p C p — p), 
that diverges at the random close packing density p cp , 
see [13]. The value for the constant b Q depends on the 
scale we choose for T g . Requiring (somewhat arbitrarily 
but adhering to convention) that ^ksT g is, in the rar- 
efied limit of granular gas, the kinetic energy of a single 
grain, we find 60 ~ 10 _25 J/(KgKelvin 2 ). • Finally, e, m 
are respectively the average internal energy and mass per 
grain. Being linear in p, W3 does not contribute to Pt 
and is not relevant for compaction. Given Eq ([3]), the 
associated pressure may be calculated from Eq ((T|) as 

P = (1 - a)BA 15 + ap 2 bT 2 g /[2{p cp - p)]. (4) 

Eqs (|2I4I) are all one needs to account for both branches 
of compaction. This ends the presentation of GSH. 

Assuming constant T gi P while evaluating Eqs (|2l4p 
simplifies the calculation and focuses the attention on 
the universal aspects of granular compaction. The un- 
derlying physics is: At given T g , the elastic deformation 
A relaxes according to Eq @ . To maintain the constant 
pressure P, the density must increase as prescribed by 
Eq (j4]). This is the irreversible branch of compaction. 
Relaxation stops at A = 0, when the pressure is purely 
temperature generated, P — Pt = ap 2 bT g /[2(p cp — p)], 
and the density assumes the final value /?/, 



2p cp /p f =! + ./! + 2ab oPcp T 2 /P, 



(5) 



obtained by taking b =const. Modifying the magnitude 
of perturbation now, pf will adjust - becoming larger 
the smaller T g is. This further density modification is of 
thermodynamic origin and reversible. 

Changing T g midway, with A still finite, leads to a 
change in A as well, if P of Eq (0| is given. This dis- 
rupts the relaxation of A, in essence resetting its initial 



condition - as observed in 14J and taken as a memory 
effect. (Generally speaking, "memory" is usually a result 
of hidden variables: When all variables have the same 
values, but the system still behaves differently, we speak 
of memory, or history-dependence. But an overlooked 
variable that has different values for the two cases will 
naturally explain the difference. The manifest and hid- 
den variables here, clearly, are p and A, respectively.) 

T g is of course not always constant, but neither does 
it have to be. It is constant for constant shear, or sound 
waves that oscillates faster than T g relaxes; it is not for 
tapping or sinusoidal shear, which give rise to periodic 
flare-ups of T g . Yet because compaction is the result of 
tiny plastic deformation accumulated over many periods, 
it is the temperature (T g ), averaged over many periods, 
that is relevant. On this longer time scale, constant T g 
is an appropriate approximation. Also, the pressure P 
is frequently nonuniform, such as when grains are filled 
in an open vessel, though it is still given. Being local 
expressions, all above equations remain valid, and we can 
read off say pf(r) from Eq ((5|) if P(r) is known. 

In a complete theory, one needs to go beyond the uni- 
versal aspects of compaction, and account for how T g is 
excited by the given perturbation and related to its am- 
plitude. We shall not do this here, only note that one can 
in principle achieve this by solving GSH for given bound- 
ary conditions. These include oscillatory shear, intense 
sound waves, periodic liquid injections, and horizontal 
tapping in various geometries, but unfortunately not ver- 
tical tapping, because particles flying ballistically part of 
the time is not a phenomenon amenable to any hydrody- 
namic description. 

The simplest perturbation to account for is oscillatory 
shear. As long as the yield surface is not breached, and 
the frequency of oscillation small compared to the relax- 
ation rate of T g , one obtains (from the balance equation 
for s g ) a quasi linear relationship between T g and the 
strain rate, T g ~ -JViftHj, see [2, 13]. Liquid injection is 
similar but involves spatial variation. A stationary sound 
field, if uniform, is again simple, as it maintains a con- 
stant T g via the strain rate that accompanies any sound. 
When the sound is turned off, T g is quickly gone. This 
increases A to maintain the pressure, and quenches the 
density at the value characteristic of T g . This is why we 
may stop the perturbation to measure the density. (Dur- 
ing the quench, there is an elastic density change that is 
negligibly small, but may be accounted for anyway.) 

If a sand-filled vessel is tapped horizontally, the blow 
generates a pulse of sound waves. As the pulse propa- 
gates through the medium and reflects multiply off the 
wall, it becomes a more or less uniform sound field - of 
brief duration but still with an accompanying T g that is 
monotonically related to the strength of the blow. And 
the quenched density measured after the tap is the one 
characteristic of this T g . All these may in principle be 
evaluated employing GSH. Tapping the vessel vertically, 
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circumstances are somewhat different, as a fraction of the 
particles will fly ballistically some of the time. During the 
flight, the grains will keep their kinetic energy, but grad- 
ually lose A: When the squeeze from the neighboring 
grains is suddenly gone after the tap, each grain oscil- 
lates around its free shape with a diminishing amplitude. 
This is puzzling, since compaction seems a succession of 
brief A-rclaxation, each picking up the work where the 
last one left it, and it is unclear why this process is not 
disrupted by A vanishing repeatedly. A conceivable an- 
swer is that these taps are not that strong. 

Newer evidence [l5| shows that the control parameter 
of the reversible branch is not the acceleration T, but 
the initial velocity tT, with t the tap duration. This is 
consistent with GSH, because with u>2 ~ s 2 g ~ T g ~ (*r) 2 , 
Eq ([5|) indeed gives pf as a function of tf. 

Eqs (|2I4[) are now evaluated with T g , P taken as con- 
stants. Starting from Eq @ and the continuity equation, 
dtp/p = —vu, we obtain (1 — a)d t p/p = d t A + A/t. 
From this, we eliminate either <9 t A or dtp by employing 
d t P(A 7 p) — 0, arriving respectively at 

d t A = -A/ta, d t p = p/r p , (6) 

with t a = t[1 + (1 - a)/pA], t p = (r/A)[l - a + pA], 
and A = (dP/dp)/(dP/dA) = - dA/dp\ P . (Employing 
Eq (U), both r p and ta may be written as a function of 
either A or p alone.) Note t p , ta are not proper relax- 
ation times, as they depend on A, p. Nevertheless, we see 
that A diminishes, while p grows - or relaxes backward 
in time. Since dP/dA ~ yAj or ta — > r for A — » 0, 
A will eventually relax exponentially toward zero, char- 
acterized by the constant relaxation rate 1/r = XT g for 
T g = constant. At the same time, p grows, ever more 
slowly, with the diminishing rate 1/t p ~ A 15 /t. It stops 
completely at p = pj for A = 0. Integrating Eqs (|6]) 
yields pit), accounting for irreversible compaction. The 
reversible branch corresponds to the change of pj with 
T g as given by Eq ([5]) . As this relation is thermodynamic 
in nature, it holds independent of path. 

We conclude that GSH is capable of a transparent, de- 
mystified account for the universal aspects of compaction. 
It describes slow density growth in the presence of gran- 
ular jiggling, with a characteristic time that diverges to- 
wards the end. This happens concurrent with the lost 
of elastic deformation. The final density, a thermody- 
namic function of T g , is larger the smaller T g is. Two 
things still need to be done: First, since a crucial start- 
ing point of GSH is the second law of thermodynamics 
using a conventional entropy, we need to understand its 
relation to the Edwards entropy. Second, we shall make 
contact with experiments, to find agreement that are, in 
spite of GSH's qualitative status, satisfactory. 

We revisit granular statistical mechanics and the Ed- 
wards entropy Secl- Substituting the volume V for the 
energy E, and compactivity X for the temperature T, 



this theory takes dV = X&Se<i as the basic thermody- 
namic relation for a "mechanically stable agglomerate of 
infinitely rigid grains at rest" [121 ]. The entropy SecL is 
obtained by counting the number of possibilities to pack- 
age grains for a given volume, equating it to e SEd . Be- 
cause a stable agglomerate is stuck in one single configu- 
ration, some perturbation is necessary to enable the sys- 
tem to explore the phase space of SecI- This ansatz, we 
believe, is best appreciated by taking the entropy more 
generally as S(E, V), or dS = (l/T)dE + (P/T)dV , with 
l/T = dS/dE, P/T ee dS/dV. For infinitely rigid grains 
at rest, we have E = 0, because the energy remains 
zero however these non-interacting grains are arranged. 
Therefore, dS = (P/T)dV, or dV = (T/P)dS = XdS. 

In GSH, grains are neither assumed infinitely rigid, 
nor always at rest. The energy, containing both elas- 
tic and kinetic contributions, does not vanish - though 
it of course does for infinitely rigid grains at rest, for 
Uij,Vi,T g ee 0. Therefore, Se<i is a minimalistic ap- 
proach, a special limit of GSH. Note, however, the fol- 
lowing objections: • Because of the Hertz-like contact 
between grains, very little material is being deformed at 
first contact, and the compressibility diverges at vanish- 
ing compression. This is a geometric fact independent 
of how rigid the bulk material is. Infinite rigidity is 
never a realistic limit for sand. • The number of pos- 
sibilities to arrange grains for a given volume concerns 
inter-granular degrees of freedom. These are vastly over- 
whelmed by the much more numerous configurations of 
the inner-granular degrees of freedom. Maximal entropy 
S implies minimal macroscopic and mesoscopic energy. 
It is quite unrelated to maximal number of possibilities 
to package grains. • Being an equilibrium consideration, 
SEd is meant to account for the reversible branch of com- 
paction. Yet Eq ([5]) is a consequence of the pressure ex- 
erted by agitated grains, not grains at rest. • Granular 
media have two different equilibria, a solid-like one for 
T g = 0, when the system is jammed, and a liquid-like 
one for finite T g , sec III of 0. A relaxing A represents 
an advance toward the liquid equilibrium, not reflected 
in SEd- For all these reasons, we are skeptical that Se<1, 
though clearly a simplification, represents a useful limit. 

When comparing experimental data to our results, 
there are some obvious difficulties: • T g is typically 
nonuniform, especially when the energy is injected from 
one side. Nonuniform shear flows, such as convection, 
will aggravate the problem. • There is boundary domi- 
nance when the system's dimension becomes a few tens 
of the particle diameter. • Crystallization may happen 
for mono-disperse particles. • The relation between T g 
and the amplitude T of the perturbation needs to be eval- 
uated for each case individually, and is usually unknown. 
Moreover, compaction experiments are executed at con- 
stant r, yet we consider compaction at constant T g . As 
T g (T) depends in general on p, A, • • ■ , the results are not 
equivalent. • Experimentally, density change is given as a 
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function of the number n of discontinuous perturbations 
such as tapping. The theory specifies the progress of com- 
paction in terms of the time t. The relation t(n) that will 
depend on T needs to be specified. Given these circum- 
stances, we shall use some of the experimental data to 
extract T g (T) and t(n), and look for agreement with the 
rest. This is admittedly much less than a complete cal- 
culation, but stays true to the spirit of accounting for the 
universal aspects of compaction. 

Denoting the initial density as po, and the final, In- 
dependent one as pf, the measured data of the Renncs 
group were given as a stretched exponential (or the 
KWW law): p f - p = (pf - p ) exp[- (n/r K ww)' 3 ], with 
/3 = 0.6 - 0.8, tkww = t x exp(r /r), r = 0.9, and 
r = 8 (relative acceleration). The Napoli group [I?} 
reported (3 = 1 for their simulation and experiments, 
and considered both results in agreement. 

On the theoretical side, the GSH parameters are the 
same we have consistently employed: a, B, bo are as in- 
dicated above. In addition, we have A = Ao (1 — p/ Pep) 
with Ao constant, and B/P = 3 x 10 6 (by taking the av- 
erage P as the hydrostatic pressure of 0.1 m and a = 0). 
Assuming an Arrhanius type law for both T g (T) and t(n), 
we also take T g (T) = er Tl/T x yj AP '/ '(2ab p cp ) (with 
A = 0.28 and T 1 = 0.56), a nd t/n = T n e~ r "/ r (with 
r„ = 7.8 and r„ = 2.2 x ltf^b oPcp /B Aq 1 ). Employing 
T g (T), Eq © becomes 2p cp /p f = 1 + Vl + Ae" 2r V r , 
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FIG. 1: Comparison between GSH and the experiments of 
[l6l. Ir?l|. (a) is the reversible branch, final density pf versus 
relative acceleration F; (b) is the irreversible ones, displaying 
X = (pf — p)/(pf — pep) versus the number n of tapping or 
liquid blows. All full lines are GSH results. The circles in (a) 
and the dotted lines in (b) are from the tapping experiment 
of [r|, while the dash lines in (b) are from liquid blowing 
experiment and simulation [l7]]. The discrepancies between 
theory, experiments and simulation are clearly of the same 
order of magnitude, see text for details. 



which fits the experimental data, see (a) of Fig [TJ Next, 
we rewrite T p (p,T g ) as T p (p,T) using T g (T) and Pf(T), 
and numerically solve the equation dtp = p/t p . The re- 
sultant p(t) —¥ p(n) for various T is shown in (b) of Fig[T] 
(The assumed Arrhanius type laws correspond to the in- 
formation contained in two of the five curves in Fig[T]) 



Summary: Compaction is a generic phenomenon, in- 
different to the method of perturbation. GSH echoes this 
by providing suitable mechanisms for its two branches, 
yielding a transparent, de-mystified account. This fur- 
ther unifies our understanding of granular behavior as 
encapsulated by GSH, which has now been shown capa- 
ble of accounting, in addition, for static stress distri- 
bution, elastic waves, elasto-plastic motion, the critical 
state, dense flow, fluidization, and jamming. 
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